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Abstract. Many, if not all, post AGB stellar systems swiftly transition from a spherical to 
a powerful aspherical pre-planetary nebula (pPNE) outflow phase before waning into a PNe. 
The pPNe outflows require engine rotational energy and a mechanism to extract this energy 
into coUimated outflows. Just radiation and rotation are insufficient but a symbiosis between 
rotation, differential rotation and large scale magnetic fields remains promising. Present obser- 
vational evidence for magnetic fields in evolved stars is suggestive of dynamically important 
magnetic fields, but both theory and observation are rife with research opportunity. I discuss 
how magnetohydrodynamic outfiows might arise in pPNe and PNe and distinguish different 
between approaches that address shaping vs. those that address both launch and shaping. Sce- 
narios involving dynamos in single stars, binary driven dynamos, or accretion engines cannot be 
ruled out. One appealing paradigm involves accretion onto the primary post- AGB white dwarf 
core from a low mass companion whose decaying accretion supply rate owers first the pPNe and 
then the lower luminosity PNe. Determining observational signatures of different MHD engines 
is a work in progress. Accretion disk theory and large scale dynamos pose many of their own 
fundamental challenges, some of which I discuss in a broader context. 
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1. Introduction 

Asymmetries were observed in planetary nebulae PNe before the Hubble Space Tele- 
scope but their ubiquity, the rapidity with which they develop, and coUimated outflow 
power have led to a cultural change in the field over the past decade (e.g. Balick & Frank 
2002). Asymmetric p/PNe have become the standard and understanding how these asym- 
metries arise is now fundamental rather than anecdotal. Perhaps most dramatic is that 
coUimated momenta of the pre-PNe (pPNe) (the reflection nebulae precursor to the ion- 
ized nebulae of PNe) exceed that which can be supplied by radiation pressure alone 
(Bujarrabal et al.(2001)). This is reminiscent of the young stellar object (YSOs) subject 



decades ago (Reipurth & Bally(2001)). The recent cultural change toward the view that 



the end states of stars (including supernovae e.g. Wang & Wheeler 2008) are asymmetric, 
offers a plethora of research opportunity. 

The beginning and ending of stars share the ingredients of in-fall, collapse, turbulence, 
and angular momentum transport. The associated increased free energy in differential 
rotation can be tapped to amplify fields and produce outflows, so common underlying 
MHD principles likely at work. Focusing on pPNc/PNc, I review the role of magnetic 
fields in these systems, and current evidence for their infiuence. I also address broader 
fundamental questions in magnetic field generation and angular momentum transport 
for which pPNe are yet another laboratory. The role of large scale fields and the need to 
connect different MHD engines with observations are central to the theme. 
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2. Basic Properties of p/PNe 

Generally, pPNe exhibit a fast bipolar outflow embedded within a slow spherically 



symmetric wind from the AGB star Bujarrabal et al.(2001) Presently, data do not rule 



out all pPNe having gone through a strongly asymmetric outflow and all PNe having gone 
through an asymmetric pNE phase. Though AGB stars produce spherically symmetric 
outflows, the initial pPNe asymmetry emerges within < 100 yr (e.g. Imai et al. 2005). PNe 
could reflect the late stages of the mechanism that produces pPNe with less asymmetry 
at very late times as supersonic motions damp. 

For pPNe ( Bujarrabal et al.(2001) I, each fast wind has a typical age At ^ 10^ — lO'^yr, 



speed ^ 50km/s, mass Mf ^ 0.5Mq, outflow rate, Mf ~ 5 x 10 ^MQ/yr, momentum 



n ~ 5 X lO'^^g.cm/s, and mechanical luminosity L,n j ^ 8 x lO'^^erg/s (can be as high as 



If ^ u.UJV'jQ, ouLuow raue, ivi j 

lO^^erg/s). The slow pPNe wind has an age Ai ~ 6 x lO'^yr, a speed ^ 20km/s a mass 
Ms ~ O.SAf©, outflow rate, Mg ^ lO^^Afo/yr, momentum lis ~ 2 x lO^^g cm/s, and 
mechanical luminosity Lm.s ^ 10^'*erg/s. For PNe, observations suggest (Balick & Frank 
2002) an age Ai ^ lO'^yr a slow wind of speed Vs ~ 30km/s of mass Ms ~ O.lAf©, outflow 
rate, Ms ~ lO~^M0/yr, momentum ~ 6 x lO'^^g cm/s, and mechanical luminosity 
Lm,s ~ 3 X 10^'^erg/s.v PNe have fast winds of speed as high as w/ ^ 2000km/s, mass 
Mf ~ 10~''A/o, outflow rate, Mf ~ lO~^M0/yr, momentum 11/ ^ 4 x 10'^''g cm/s, and 
mechanical luminosity Lmj ^ 1.3 x lO'^^erg/s. 

The pPNe phase demands the most power and the linear momenta of fast bipolar pPNe 



outflows seems too large for radiation driving Bujarrabal et al.(2001) This motivates the 



need to tap rotational energy, either from redistribution of angular momentum within a 
single star or via deposition of angular rotational energy from binaries. 



3. Observations of Magnetic Fields 

Detection of magnetic fields in p/PNe is not by itself a proof of their dynamical impor- 
tance as the relative strength of the field and the local kinetic energies must be considered. 
Complementarily, a weak magnetic field at large distances from the outflow engine does 
not rule out a magnetically dominated launch at the engine. 

Magnetic field detection is technically challenging but careful analysis has led to esti- 
mate of magnetic field strengths and geometries in several important classes of late AGB 
and post- AGB systems by different techniques. In water maser sources such as W43 (an 
evolved AGB star with a precessing jet signaturing the beginning of the pPNe stage) 
magnetic fields of 85mG have been inferred at radii of 500AU from the engine via in- 
terpretation of circular polarization of H2O masers (Vlemmings et al. 2006). The 85mG 
fields at 500AU and their inferred geometry seem consistent with what is needed to dy- 
namically influence the flow. (See Sabin et al. 2007 for field geometry in other sources). 
The precession suggests a binary interaction, but whether the field is dynamo produced 
and/or connected to a companion within the envelope or the AGB core is uncertain. 
Vlemmings (2007) reviews maser field measurements of giant stars highlighting that the 
H2O masers probe field on ^ 10^ AU scales, OH masers probe 10'^ AU scales and SiO 
masers probe ^ lOAU. The complied measurements show that, statistically, fields fall off 
faster than but slower than r^^. This is a weak constraint on models. 

From VLT spectropolarimeteric analysis of several central stars of PNe (CSPN) , Jordan 
et al. (2005) detected magnetic field strengths of order kG. The scale of the photosphere 
for these system is ^ 2 x 10^°cm. These field strengths are are consistent with the larger 
scale scalings presented in Vlemmings et al. (2007) and the strengths of fields needed to 
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power the pPNe outflows by Poynting flux, and fields produced by dynamos in the AGB 
engines (Nordhaus et al. 2007). 

Soker and Zoabi (2002) and Soker (2006b) have argued for weaker fields and rather 
than a primary driver of outflows, an indirect shaper of outflows possibly via influence on 
the dust distribution and geometry of radiative driven winds. But this would leave un- 
solved the large coUimated momenta of pPNe engines Bujarrabal et al.(2001) ) for which 
radiation driving is insufficient. 



4. Dynamical Magnetic Shaping and/or Launching 

The "launch" region of MHD outflows is the region where the magnetic force dominates 
the flow. This extends to a height typically no greater than Zc 50Ri, where Ri is the 
inner-most radial scale of the engine, (e.g. the inner radius of an accretion disk). In the 
launch region the bulk flow is sub-Alfenic below Zc- The "propagation" region describes 
z > Zc where the poloidal flow speed exceeds the Alfen speed, eventually approaching its 
asymptotic speed. Presently, only the propagation regions are observationally spatially 
resolved. Also, because there are ~ 2 or more orders of magnitude between the engine 
where dynamos and jets formation operates and the asymptotic propagation region the 
computational demands prohibit simulatating the combined physics of the launch region, 
jet formation and asymptotic propagation. As emphasized in sec 7., even the nonlinear 
physics of just an accretion disk has been prohibitively computationally expensive. We 
must patch together different pieces of physics from different approaches and scales to 
extract a complete picture. 

4.1. Shaping with Imposed Fields 

In a semi-analytic approach. Chevalier and Luo (1994), imposed a magnetic field geom- 
etry such that the toroidal field falls off more slowly than the poloidal field and at large 
distances and a spherical wind shocks the ambient gas increasing the toroidal field via 
compression. For their parameters, the magnetic shaping occurs only after the shock. 

Garcia-Segura (1997) and Garcfa-Segura et al. (1999) were the first to simulate the 
direct effect of imposed toroidal fields on PNe shaping by driving a hydrodynamic wind 
into a pre-magnetized medium. The field strengths were consistent with those of sec. 
3 above, and their infiuence (along with stellar rotation, and photoionization) on the 
shaping was summarized in simple schematic (Garcfa-Segura 1999) . Rotation -|- magnetic 
fields leads to strong coUimation compared to the weaker coUimation of just the rotational 
influence on radiatively driven winds. Gardiner & Frank (2001) discussed that the shaping 
observed in such simulations were outside the restrictive parameter regime of Chevalier 
& Luo (1994) and that the influence of reasonable strength fields should generically affect 
the outflow before the shock compression. . 

4.2. Launching and Shaping with Imposed Initial Fields 

Because the observations require both a mechanism of launch and shaping, it is neces- 
sary to consider the role of magnetic flelds in both. Garcfa-Segura et al. (2005) imposed 
a toroidal fleld the surface of the AGB star to launch an outflow via the fleld pressure 
gradient and studied the propagation and asymptotic coUimation over 2 orders of magni- 
tude in scale. Given field strengths of ^ 40G (consistent with observations) at the surface 
of the AGB star, outflows consistent with the required power and shaping were observed. 

Matt et al. (2006) added additional physics in demonstrating how an initially weak 
field can grow linearly in time by extracting rotational energy of its anchoring base, and 
produce a rapid bipolar outflow from the gradient in toroidal field pressure. They took 
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a gravitating spherical core surrounded by an initially hydrostatic envelope of ionized 
gas. A dipolar magnetic field was anchored on the core, threading the envelope. The core 
was set to rotate initially at 10% of the escape speed. The toroidal field amplified from 
the differential rotation between core and envelope and when the toroidal field pressure 
gradient overcame the envelope binding energy, envelope material was rapidly expelled 
in a quadrupolar outflow: The wound up diploe field had maximum toroidal field at 
intermediate poloidal angles in each hemisphere, so material was both along the poles 
also squeezed out from the equator. 

4.3. Launching via Dynamo Engines 

The models of the previous two subsections do not fully address where the dynamically 
important large scale fields come from. Even the Matt et al. (2006) paper invoked linear, 
laminar field growth and treating the rotator base as a boundary condition. 

Papers such as those by Pascoli (1993,1997) Blackman et al. (2001a, b), Nordhaus et 
al. (2007) use flavors of a mean field dynamo theory for stellar or disk engines to estimate 
field strengths and Poynting fluxes that could arise. These papers are not simulations and 
do not track global field geometry, or follow the production and evolution of outfiows 
from the field production region. Furthermore, the nonlinear physics of dynamos (see 
sec. 7) require approximations when cmedding them in astrophysical scenarios. Despite 
their shortcomings, such toy models do produce promising results with respect to field 
strengths and Poynting fiuxes 



5. Single Star vs. Binary Models 

5.1. Single star models not yet ruled out 

Binaries can easily supply the needed free energy to amplify fields via differential rotation 
(ordhaus & Blackman 2006). Isolated star MHD models cannot be ruled out, but there 
are caveats (Soker 2006b; Nordhaus et al. 2007) as I now discuss. 

Blackman et al. (2001) investigated an isolated interface dynamo model operating 
at the base of the AGB convective zone in which angular momentum is conserved on 
spherical shells as the star evolves off the main sequence and the resulting rotation profile 
provides the available differential rotational energy from which the field is amplified. To 
drive bipolar pPN/PN, the corresponding dynamo must operate through the lifetime of 
the AGB phase (lO^yr) until radiation pressure has bled most of the envelope material 
away. Only then can the Poynting flux unbind the remaining material. But the dynamo 
would drain differential rotation on time-scales short compared to the AGB lifetime. Only 
if differential rotation is extracted from deep within the core or re-seeded by convection 
can the drain be overcome (Nordhaus et al. 2007). Viability of the single star dynamo 
outflow model, depends on whether or not the differential rotation can be re-supplied 
over the AGB lifetime via convection as in the Sun. 

5.2. Common envelope binary scenarios 



In a common envelope (CE) Iben & Livio (1993) the companion drags on the envelope 



of the primary, transferring angular momentum and energy. If the envelope cooling time 
is long enough, a fraction a ^ 0.1 of the loss in gravitational energy can spin up and 
unbind the envelope. |Reyes-Ruiz fc Lopez (1999) discuss the secondaries for which ac- 



cretion disks will form around the primary core. Brown dwarf (BD) (O.OOSM© < Merit ~ 
O.O7M0) radii increase with decreasing mass while their Roche radii decrease with de- 
creasing mass. Such objects unstably lose mass. Since the circularization radius lies 
outside of the primary's core, a disk can form within a few orbits. This contrasts the 
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M > Merit case for which the stehar radius decreases with decreasing mass more strongly 
than the Roche radius. Supercritical companions have a circularization radius within the 
primary's core. Material leaving the secondary would then initially spiral swiftly into the 
primary rather than orbit quasi-stably. Nordhaus & Blackman (2006) show that accretion 
from planets and low mass stars may also be important (and perhaps more common be- 
cause of the brown dwarf desert Grether & Lineweaver(2006) ) and can supply the needed 
accretion power for pPNe. Planets are of sub-critical mass and a disk could form as per 
a BD, before unbinding the envelope. 

A companion star is supercritical so if stellar incurs Roche lobe overflow after CE, 
sustained accretion requires the binary to lose angular momentum. This could happen 
as the overflowing secondary drags on residual inner envelope material. Even though 
the circularization radius for M2 > Merit is inside the core, the in-spiraling material still 
incurs differential rotation and could amplify magnetic fields. A significant energy release 
via accretion in the AI2 > Merit case could occur on impact to the inner core, producing 
dwarf novae bursts masked by the stellar envelope. If this initial accretion phases can drop 
M2 below Merit, keep M2 filling its Roche lobe, and leave enough angular momentum to 
form a Keplerian disk, then accretion could proceed as for the initial M2 < Merit case. 



Huggins (2007) suggests a ~ O(IOO) year delay between ejection of cicumbinary dust 



tori and jets in pPNe. If a fraction of ejected envelope material becomes the dust torus, the 
delay could be the time for the companion to lose enough mass to move the circularization 
radius outside the core for the supercritical case, or a viscous time. 



6. Accretion Disk Outflows in pPNe and PNe: 



Accretion disk outflows have a mechanical power Blackman et al. (2001b) 



GM.Mge 36 M,M_4 
=4.5 X lO^^e-i— , (6.1) 

where e is the efficiency of conversion from accretion to outflow, Ri is the inner disk 
radius. Ma is the accretion rate, G is Newton's constant, M* is the central stellar mass, 
i?j,io = R^/10^°cm, e_i = e/0.1 and Ma -4 = Ma / 10''^ Mq/jt. For an MHD outflow, 
Eq. (|6.ip equals the Poynting flux at the launch surface. Propagation into a low density 
medium produces an asymptotic outflow speed ^ Q,rA where is the angular speed 
of field anchor point and is the radius where the poloidal outflow speed equals the 
Alfven speed. This product is typically a few times the escape speed of the inner most 

radius of a disk and is thus at least Vout Vesc = 1600 ( jj^'^ ) km/s 

6.1. Accretion onto Primary 

The Vout above depends only weakly on Ma via Ri, but strongly on the inertia of material 
blocking the outflow: Momentum conservation gives 

jaMenv + Mf 

where Mf/M^ — 3.3 x 10~''e_iMao,-3 /j^^"" r^^/^dr is the mass in one of the fast col- 
limated outflows, Menv is the envelope mass, /o ~ 0.2 is the solid angle fraction inter- 
cepted by the collimated outflow, and r = t/lyr is used to incorporate Ma oc t^^^^ of 



Reyes- Ruiz & Lopez (1999) The numbers are scaled to pPNe so Menv » Mf and for 
an envelope of mass 2Mq, the intercepted mass is 0.2 Mq for /o = 0.2. 

Eq. (|6.2p is the observed speed of the fast when blocked and loaded by the envelope. By 
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the end of the pPNe phase, the envelope is quite extended, reducing the optical depth and 
revealing material moving at the "free streaming" fast wind speed. Assuming a dust-to- 
gas mass ratio of 1/100 and micron sized grains of density of 2g/cm'^, the optical depth 

from dust is r, ~ 2.5 x lO'^ (^^-^^li^^) (to^) (to^) ' scaled for PNe. For 

pPNe, the density increases by ^ 10'* and R is down by a factor of 10, so Td ^ 2.5. The 
different optical depths of pPNe and PNe can thus explain why observed PNe fast winds 
can have Vf = Vout > 1600km/s, whilst those of pPNe have Vf — Vout < lOOkm/s. 

Time dependent accretion outflows described with Eqs. (|6.ip and (|6.2p are consistent 
with the high pPNe outflow mechanical luminosity and the fast PNe wind speed of Sec. 1 



when M^,/R^ corresponds to a WD. Reyes-Ruiz & Lopez (1999) considers a companion 



of mass M2 ~ O.OSMq < M^rit and a Shakura-Sunyaev Shakura & Syunyaev(1973) 
viscosity parameter Uss ^ 0.01, for which the accretion rate then decays as Ma ~ 1.6 x 
lO-^t-s/^Mg/yr.. Using this in with e = 0.1 for t = 100 yr with i?^ = 2 x lO^cm 
and Ml = O.GM© gives L„j ~ 4.3 x IQ'^'^ {t/lyi)-^/^ . This provides the needed power 
demands of Sec. 1 for pPNe after 1000 yr and for PNe after 10** yr. 

A more careful analysis of the predicted jet speed evolution is needed for specific outflow 
models as the envelope evolves in order to test the idea contained within the rough 
estimates just discussed. See also Garicia-Diaz et al. (2008) for a non-accretion based 
time dependent comparison of jet outflow speed from simulations with observations. 

6.2. Accretion onto Secondary 

It is also possible for accretion disks to form around the secondary (Soker & Livio 
1994; Mastrodemos & Morris 1998; Soker 2005). The Bondi wind accretion rate is 

~ (^) [i+tviv\'^Y/^ ' '^^^'^6 V is the orbital speed of the secondary and Vyj is 
the slow wind speed from the primary. In general, for M2 < Mi, reasonable parameters 
provide an accretion rate compatible PNe luminosities if the companion is either main 
sequence or compact star, but the ubiquity of high coUimated fast wind pPNe powers 



( |Bujarrabal et al.(2001) ) would require an overabundance of accreting WD companions. 



7. Key Issues in Large Scale MHD Dynamo and Accretion Theory 

7.1. Large Scale Dynamo Theory: Recent Developments, Open Questions 

Large scale dynamo (LSD) theory describes the sustenance of magnetic fields on time / spatial 
scales larger than turbulent scales. Whether large scale fields for jets are advected or LSD 
produced has been debated but but ultimately, the equations that include a the com- 
petition between turbulent transport, fiux accretion, as well as LSD action need to be 
solved. Field reversals in the sun prove that an LSD can and must operate therein. 

For ^ 50 years, a problem with textbook LSD theory (e.g. Moffatt (1978) ) has been 
the absence of a proper saturation theory that predicts how strong the large scale fields 
get before non-linearly quenching via the backreaction of the field on the driving flow 
kicks in. But substantial progress toward a nonlinear mean fleld theory has emerged 
in the last decade via a symbiosis between analytical and numerical work. Coupling 
the dynamical evolution of magnetic helicity into the dynamo equations turns out to 
be fundamental for predicting the saturation seen in simulations. For recent extensive 
reviews see (Brandenburg & Subramanian 2005; Blackman 2007). 

Much of the work in LSD theory, has focused on systems that are initially globally 
reflection asymmetric (GRA). This means a global pseudoscalar is imposed by the bound- 
ary conditions-for example, rotation and stratification lead to the kinetic helicity pseu- 
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doscalar, common to the standard textbook (Moffatt (1978)) "adyn effect" of mean field 



dynamos. However magnetic helicity is actually the unifying quantity for LSD. There 
are two classes of GRA LSD for which an electromotive force aligned with the mean 
magnetic field is essential (see e.g. Blackman 2007). The first is fiow driven helical dy- 
namos (FDHD) which occur inside of astrophysical rotators. Here the initial energy is 
dominated by flows and the field responds. These are linked to a corona by magnetic 
buoyancy. In coronae, the second type of LSD, the magnetically driven helical dynamo 
(MDHD) can operate. This characterizes relaxation of magnetic structures to larger (jet 
mediating) scales in a magnetically dominated environment subject to the injection of 
smaller scale magnetic helicity. The MDHD is directly analogous to laboratory plasma 
dynamos that occur in reverse field pinches (RFPs) and Spheromaks. 

LSDs always involve some helical growth of the large scale field which is coupled to 
a helical scale fields of opposite sign. When small scale magnetic or currently helicity 
evolution is coupled to the large scale field growth, the modern mean field 'dynamo adyn 
effect that predicts the correct satuation becomes the difference between kinetic helicity 
and current helicity: For a FDHD simulated in a closed box (Brandenburg 2001), the 
current helicity builds up as the large scale field grows and quenches the FDHD (Black- 
man & Field 2002). Complementarily, for an MDHD, the system is first dominated by 
the current helicity and a growing kinetic helicity then acts as the backreaction (Black- 
man & Field 2004). Both FDHD and MDHD are accessible with in the same framework, 
all unified by tracking magnetic helicity evolution, and aided by thinking of the field as 
ribbons rather than lines (Blackman & Brandenburg 2003). More work on how the fields 
evolve from within the rotator to produce the global scale fields in coronae which in turn 
produce jets are needed as most work on dynamo theory has focused on the FDHD. 

Because the buildup of small scale magnetic (or current) helicity quenches the LSD, 
preferential ejection of small scale helicity vs. large scale helicity through a boundary 
can in principle alleviate this quenching (Blackman & Field 2000; Vishniac Cho 2001; 
Sur et al. 2007). Numerical simulations support this general notion, particularly when 
shear is present and when surfaces of shear align toward open boundaries (Brandenburg 
& Sandin 2004; Kapyla et al. 2008), thereby allowing needed helicity fluxes. 

LSD action has also been observed in non GRA simulations (e.g. Yousef et al 08; Lesur 
& Ogilvie 08) implying that the minimum global ingredients for this class of analytic 
LSD is shear, plus turbulence that feeds azimuthal held back to toroidal fleld. There is 
work in progress to understand the simulations guided by analytic models (e.g. Vishniac 
& Brandenburg 1997; Blackman 1998; Kleeorin & Rogachevskii 2003; Schekochihin et al 
2008). The non GRA LSDs grow large scale flelds on scales larger than the turbulence but 
smaller than the global scale. In coherence regions, there is a field aligned electromotive 
force (EMF), and thus an intermediate scale source of magnetic helicity that may switch 
signs between coherence regions and globally average to zero. It may be that the non- 
GRA LSD action always involves a local helicity flux between coherence regions. 

7.2. Accretion Disks: More questions and the need for large scale fields 

The magneto-rotational instability (MRI) has emerged as a leading candidate for angular 
momentum transport in accretion disks (e.g. Balbus & Hawley 1998) Although the MRI 
exists without the LSD, they are likely coupled in nature. I explain this below. 

There is a disconnect between what shearing box simulations have told us about the 
MRI vs. how the instability might operate in nature. To date, simulations have primarily 
told us that the MRI is plausible but do not produce a robust theory of saturation or 
robust values of transport coefficients for modelers. This may frustrate, but patience 
(for several more decades) is required, as the computational and conceptual demands 
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are substantial. For example, to achieve better resolution, most first generation MRI 
simulations (except Brandenburg et al.l995) did not use explicit viscosity or magnetic 
difFusivity (see discussion in e.g. Fromang et al. 2007). the magnetic Prandtl number 
affects the magnetic spectrum and the transport coefficients. In adddition, the value 
of the angular momentum transport coefficient ass depends strongly on the box size 
and the strength of the initially imposed weak mean field strength (Pessah et al. 2007). 
Interestingly, ass varies ~ 4 orders between simulations, but assP, where f3 is the ratio 
of thermal to magnetic pressure, is nearly a constant (Blackman et al. 2007). 

In addition to the need for explicit diffusivities, perhaps the most important fron- 
tier is actually role of large scale magnetic fields in angular momentum transport and 
thus, non-local contributions to magnetic stresses that transport angular momentum. Its 
importance is motivated from three different paths: 

(1) Large scale fields are evidenced from theory, simulation, and observations of coronae 
and jets. Large field structures more easily survive the buoyant rise to coronae without 
being shredded by turbulence within the disk. Plausibly, the integrated stress for struc- 
tures on scales above that which survives the buoyant rise would contribute to the large 
scale magnetic stress. But shearing boxes that study only a local region non-local large 
scale magnetic stress is excluded. In a real system this could be the most important part. 

(2) Shearing box simulations artificially impose a steady-state because differential rota- 
tion is imposed as a steady forcing not subject to the backreaction of the amplified field. 
Hubbard and Blackman (2008b) , argue that this may be more restrictive that previously 
recognized: In a real disk, energy in differential rotation is susatined only by accretion 
itself. If a steady-state is to be maintained via turbulent transport alone then there must 
be 100% power throughput from differential rotation to the turbulent cascade. This is 
not guaranteed if large scale fields drain power. An accessible steady-state solution must 
then incorporate stresses from large-scale fields in addition to turbulent transport. 

(3) The role of large scale magnetic fields is consistent with the implications of (Pessah 
et al. 2007) which shows that MRI stresses in simulation boxes where the radial extent 
is of order the vertical scale height scale with the ratio of box size to scale height. The 
contribution of large scale fields would increase as the radial and vertical scales are in- 
creased, highlighting a strongly non-local contribution to stresses that transport angular 
momentum. Evidence for non-local MRI behavior is also seen in Bodo et al.(2008) Box 
sizes for thin disks must be extended in the radial direction and accordingly in the vertical 
direction as buoyant loops tend to have radial scales comparable to vertical scales. 

Ultimately, mean field accretion disk theory should be coupled to an LSD theory in a 
real disk as they are actually artificially separated components of what should be a sin- 
gle mean field theory. Note that for any turbulent disk, any assumption of axisymmetry 
for bulk quantities automatically requires that theory to be a mean field theory. This is 
often veiled by a mere replacement of the actual viscosity with a turbulent viscosity as 
in the Shakura-Sunyaev approach (Shakura & Syunyaev(1973) I. Hubbard & Blackman 
(2008a) for example, suggest that in a formal mean field theory, a term involving fluctu- 
ations of density and velocity might be interpreted as an additional transport coefficient 
in the mean surface density equation, restricting available steady states. Balbus et al. 
(1994)incorporated this term into a redefinition of accretion rate. 



8. Much work needed to link theory and observation 

It remains a major challenge to rigorously couple the engine physics of field generation 
and accretion to jet launch and jet propagation in a unified theory or simulation that 
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make distinct observational predictions for specific engines. This enterprise spans several 
subfields of theoretical astrophysics, let alone the specific application to p/PNe. 

Here however, is a list of possible lower hanging fruit for linking theory with observa- 
tions:(l) Evaluate the kinematic constraints/predictions of outflows from the scenario of 
accretion onto the primary and compare the distribution of inferred fast outflow speeds 
to what would be expected from known binary statistics of low mass stars. This would 
help constrain the commonality of accretion onto the primary vs. secondary as the lat- 
ter has a broader range of masses. (2) The more massive the companion, the more CE 
models would predict mostly Oxygen rich rather than Carbon rich post AGB systems 
because the binding energy for the early AGB phases is higher. Low mass companions 
like planets may terminate the AGB only in the thermal pulse phase. (3) Crystalline dust 
in post- AGB systems can be produced if a binary induced spiral shock anneals silicates 
?. Is this universal? (4) CE evolution would predict equatorial outflows that precedes 
any accretion driven poloidal jet. Is this consistent with the delay of Huggins (2007) and 
the geometry of equatorial outflows? (6) Are fast outflows contaminated by material that 
could represent accretion disk residue of shredded low mass companions? (7) Are time 
scales of observed outflow precession consistent with the gravitational influence of a bi- 
nary? (8) Can double peaked line profiles be detected to identify accretion disks within 
the launch region? (9) Can shrouded novae outbursts from a M2 > Merit companion 
feeding the primary be detected in X-rays? (10) Improved statistics on the fraction of 
bipolar pPNe, the fraction of suitable precursor binaries for CE, and the fraction of stars 
which evolve to be pPNe will improve evaluation as to whether all PNe incur asymmetric 
phase. (12) Can the approach of Ferreira et al. (2006) be generalized to pPNe outfiows? 
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